hep-ph/0505168 

Retardation effects in the rotating string model 

Fabien BuisseretQ and Claude Semajll 

Groupe de Physique Nucleaire Theorique, Universite de Mons-Hainaut, 
Academie universitaire Wallonie-Bruxelles , 
Place du Pare 20, BE-7000 Mons, Belgium 
(Dated: February 2, 2008) 

}Q Abstract 
O 



o 

(N 



A new method to study the retardation effects in mesons is presented. Inspired from the covariant 



oscillator quark model, it is applied to the rotating string model in which a non zero value is allowed 
for the relative time between the quark and the antiquark. This approach leads to a retardation 
term which behaves as a perturbation of the meson mass operator. It is shown that this term 
. preserves the Regge trajectories for light mesons, and that a satisfactory agreement with the 

experimental data can be obtained if the quark self-energy contribution is added. The consequences 

vo . 

of the retardation on the Coulomb interaction and the wave function are also analyzed. 
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I. INTRODUCTION 



The retardation effect between two interacting particles is a relativistic phenomenon, due 
to the finiteness of the interaction speed. The light mesons are typical systems in which 
the retardation mechanisms can significantly contribute to the dynamics, since the light 
quark can move at a speed close to the speed of light. To take into account this effect 
into effective models, one has to work with a fully covariant theory. The most elegant 
approach, from a theoretical point of view, is the formalism with constraints □ D , but it 
quickly leads to complex equations, uneasy to deal with if one wants to get analytical or 
numerical results. Another covariant approach of mesons is the covariant oscillator quark 
model (COQM), which allows to find an analytical expression for the wave function and 
numerical results in good agreement with the data 3,0]. Unfortunately, this model uses a 
quadratic potential, this form being different from the linear potential, commonly assumed to 
describe the confining part of the interaction and suggested by lattice calculations. Moreover, 
vanishing quark masses are not allowed in this approach. 

Apart from these two approaches, most of the effective models found in the literature are 
based on the equal time ansatz, which simply takes the time coordinates of both particles 
to be equal, neglecting the retardation effects. This procedure allows to deal with simpler 
equations, and relativistic corrections can be obtained by developing an expansion in v 2 /c 2 
of the model considered, like the Bethe-Salpeter equation or even the QCD Lagrangian 0,1^. 

The model we propose here is an attempt to include retardation effects into the rotating 
string model (RSM) [3, 0] without making such an expansion, in order to estimate the 
retardation contribution for light quarks. The RSM is an effective model derived from the 
QCD Lagrangian, describing a meson by two spinless quarks linked by a straight string. It 
las been shown that the RSM was classically equivalent to the relativistic flux tube model 



^, This last model, firstly presented in Ref. [ill 



121 ]. yields meson spectra in good 



agreement with the experimental data 13]. Our method, inspired from the COQM, relies 
on the hypothesis that the relative time between the quark and the antiquark must have a 
non zero value. In our framework, the evolution parameter of the system is not the common 
proper time of the quarks and the string, but the time coordinate of the center of mass 
which plays the role of an "average" time. 

Our paper is organized as follows. In Sec. |TIJ we present the general formalism of our 
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approach. We compute the retardation contribution to the meson mass in Sec. Mil We find 
different approximations for this contribution in Sec. IIVI with a special interest for light 
quarks, and we numerically study the retardation effects in Sec. El As our model relies on 
unusual hypothesis, it is worth comparing our results with those of other existing models. 
This is done in Sec. IVII Finally, we compare the meson spectra of our RSM including the 
retardation term with the experimental data in Sec. lVIII Some concluding remarks are given 
in Sec. IVII II The appendix contains some useful formulas. 



II. THE ROTATING STRING MODEL WITH NON ZERO RELATIVE TIME 

It has been shown that, starting from the QCD Lagrangian and neglecting the spin 
contribution of the quark and the antiquark, the Lagrange function of a meson can be built 
from a Nambu-Goto action 7] which reads (77 = diag(H ) and h = c = 1) 

C(t) = — mi\Jti{ — m 2 \jx^\ — a J d(3\J (ww') 2 — w 2 w' 2 . (1) 

The two first terms are the kinetic energy operators of the quark and the antiquark, whose 
current masses are m\ and 7722. These two particles are attached by a string with a tension 
a. Xi is the coordinate of the quark i and w is the coordinate of the string, w depends on 
two variables defined on the string worldsheet: One is spacelike, /3, and the other timelike, 
t. Derivatives are denoted w' = dpw and w = d T w. In this picture, r is a common proper 
time for the string and the quarks. Introducing auxiliary fields to get rid of the square root 
in the Lagrangian (P) and making the straight line ansatz to describe the string, an effective 
Lagrangian can be derived jq] 
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mi m 



— + — + a x R 2 + 2a 2 Rr - 2(c x + (a^Rr 



-2(c 2 + (a 2 )rr + a 3 r 2 + (a 4 + 2( Cl + ^a^r 2 



(2) 



where the coefficients ai, a 2 , . . . , given in the appendix, depend on various auxiliary fields 
/ii, fi 2 , v, and m The parameter ( defines the position R of the center of mass: R = 
(x\ + (1 — ()x 2 . r is the relative coordinate r = X\ — x 2 . The auxiliary fields /ii and \i 2 
are seen as constituent masses for the quarks, and v can be interpreted in the same way as 
an effective energy for the string whose "static" energy is ar P. Iic|. Let us note that the 
straight line ansatz for the string implies w = R + {(3 — Qr. 



The total and relative momentum, computed from the Lagrangian (J2J), are respectively 

dC 

PfM = — h— = -axRp - a 2 r fl + (ci + C«i)^, (3a) 
dC 

Pti = ^7 = ~a 2 R l i - astft + (c 2 + C^V (3b) 

As we will work in the center of mass frame, the total vector momentum P of the system 
must vanish, which implies that 

^ = (ci + Cai)r -a 2 r 

Moreover, the relative vector momentum p is given by 

p = a 2 R + a 3 r - (c 2 + (a 2 )r. (5) 

Thus, we impose a 2 = in order that p does not depend on the motion of the center of 
mass. This leads to the following value for the parameter ( 

^1 + ^2 + Jq d/3 v 

Equation © reduces to ( — 1/2 in the symmetrical case (mi = m 2 ) and to £ = m 1 /(mi+m 2 ) 
in the nonrelativistic limit, as expected jlfl| . 

To go a step further, one usually takes the temporal coordinates of the quarks and the 
string to be equal to the common proper time r, this time being also the time t in the center 
of mass frame 

x \ = xl = W° = T = t. (7) 

Then, we have r — (0,f), R — (t, R), r — (0,r), and R = (1,R). This procedure allows to 
deal with simpler equations, but neglects the relativistic retardation effects due to a possible 
non zero value of the relative time r°. Since these effects are precisely those we want to 
study in thispaper, we have to make a less restrictive hypothesis. As in the formalism of 
the COQM |3j,y, we define 

r = (<7,f), R=(i,R). (8) 

The temporal coordinate of the center of mass, i, can be seen as an "average time" for the 
meson. This is particularly clear in the symmetrical case, where i = (xi + x 2 )/2. Our choice 
is to take t as the evolution parameter of the system. We identify it as the common proper 
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time for the quarks and the string, and the dotted quantities are derived with respect to t. 
We have for example 

r = (&,?), R=(1,R). (9) 

The special case a = is equivalent to the relation 0. 

The Lagrangian (j2J) can now be rewritten using formulas (jlj), (JHJ), (jHJ), and (jHJ) as 



C = C + AC, 



(10) 



with 



Co 
AC 



TYli-i TYlir) X jf . Q . ; Q * > * ; 

1 h ai H (c x — a 4 ai)r — a 3 air + 2aiC 2 rr 

/ii /i 2 ax V 



(ci + C^ijcr + c 2 ao- - —b - -(a 4 + 2(,ci + C aijcr 



(11) 

(12) 



We have gathered the relative time dependent terms in AC, which contains the contribution 
of the retardation. Let us remark that the Lagrangian Co does not depend on (. 

We shall consider in the following AC as a perturbation of Co- With this hypothesis, the 
auxiliary fields can be eliminated by considering only the constraint 5Co = 0. In Ref. flo| . 
it is shown how a set of three equations can be derived from the Lagrangian ([lljh to define 
the usual rotating string model (RSM) 
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with 



F( yi 



H (arcsmyx + arcsint/ 2 J, 

Vt 



arcsinyi - y { \ 1 - y. 



and y t = yi + y2- 



(13a) 
(13b) 

(13c) 
(13d) 



p r is the radial momentum and yi can be seen as the transverse velocities of the quark i. 
The first relation gives the cancellation of the total momentum in the center of mass frame, 
while the two last ones define respectively the angular momentum and the Hamiltonian. As 
we can see in Eqs. (|13|). the only remaining auxiliary fields are //j. The extremal values of 
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the auxiliary fields rj and v are given by 

m = k([3-<J)), (14a) 
ar , , . 

v 1 - vt\P - 

with k = — <b = —^r — , and a = m , M2 . Let us remark that a closed form can not be 
obtained for the Hamiltonian H Q because it impossible to eliminate analytically the variables 
?/i and y 2 as a function of L, by means of the two first Eqs. (fT3|) . 

Since we have a contribution from the relative time, the total Hamiltonian is given by 

H = H + AH, (15) 

with 

AH = T,a - AC. (16) 

With £ = dC/da, Eq. © leads to 

a = — , (17) 

a 3 

and finally we obtain 

AH = + ^Sa - ( Cl + Cai)a - ^U 2 + \{a A + 2Cd + CV)^, (18) 

the perturbation of the RSM Hamiltonian due to the retardation effect. 

In the following, for simplicity, we will focus on the symmetrical case. Then ( = 0, and 
the Hamiltonian (fTB*|) becomes 

AH = — — [£ 2 - 2c 2 £(x + (c 2 - a A a 3 )a 2 ] . (19) 
2a 3 

Wether AH is really a perturbation or not has to be verified. We will check this hypothesis 
a posteriori by a numerical computation of the retardation contribution to the meson masses 
(see Sec. El). 

III. RETARDATION EFFECTS 
A. Contribution to the mass 



Up to now, we were working in a classical framework. But in order to study the influence 
of the retardation on the meson spectrum, we have to consider a quantized version of the 
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total RSM Hamiltonian H + AH. We can thus replace L by \J l{l + 1) and consider r and 
a as operators such that 

[r, p r ] = [a, E] = -z. (20) 
The total Hamiltonian has schematically the following structure 

H(a,r) = H Q {r)+AH(a,r). (21) 

The relative time a only appears in the perturbation, and Hq only depends on the radius r. 
So, we make the following ansatz to write the total wave function 

Wr)> = \R(f)) ® |A(a)> , (22) 

where |i?(f )) is a solution of the eigenequation 

H (r) \R(f)) = M \R(f)) . (23) 



[14]. 



Such a problem can be solved, for instance, by the Lagrange mesh technique |14| . The total 
mass is written 

M = M + (A(a) | ® (R(f) \ AH{r, a) \R(r)) ® \A(a)) 

= M + AM. (24) 

the contribution AM is then given by the solution of the eigenequation 

AH(a) \A(a)) = AM \A(a)) , (25) 

where 

AH(a) = (R(f) | AH(r, a) \R(f)) . (26) 
In order to eliminate the unphysical excitations of the relative time, we consider only the 



ground state of the Hamiltonian A7i(a), as it is done in Refs. 
this Hamiltonian is defined by 



Using formula (fT§j). 



where all mean values () are computed with a space function R(f) and where {A, B} = 
AB + BA. We also use the approximation (1/x) ~ !/(#)■ 
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Using Eqs. (jl4a|) and (jA.3|) . we see that C2 oc k. Since we are in the symmetrical case, we 
can assume («) = 0, and so (c 2 ) = 0. On the contrary, (c|) 7^ because (p 2 ) > 0. Finally 
the Hamiltonian A7i takes its final form 

AH = [S 2 + (c 2 - a 4 a 3 )(7 2 ] . (28) 

With our approximations, the retardation contribution to the Hamiltonian looks like an 
harmonic oscillator for the canonical variables (er, £) . Let us define 

p 2 = c 2 , — a 4 a 3 . (29) 

Thanks to formulas (JA.lj) . we can compute p 2 

p2 = 2n/ (2 + ^ ("^V 71 - + arcsin yj^j (y^Jl-y 2 + arcsin 

2 

" 4^ (-2/V / l-2/ 2 + arcsin y) . (30) 



Assuming that (p 2 ) > (this is checked in Sec. II V|) . the ground state solution of the 
eigenequation is given by 

AM = --a;, (31a) 

Mr) = (^*v(AA> ( 31b ) 

with 



7T/ V 2 



/? = v^>, (31c) 

w - ^ (aid) 

An immediate conclusion to draw from Eqs. (|31|) is that the retardation effects bring a 
negative contribution to the meson masses, and that the more probable value for the relative 
time is a = 0. It is worth noting that these results are formally identical to those of the 
COQM. In Sec. IVI[ we compare the two approaches with more details. 



B. Modification of the coulomb term 



In the RSM, the string contribution takes into account the interactions at large distances, 
which are responsible for the confinement. To make more realistic models, it is necessary 
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to add short range potentials For instance, the one gluon exchange mechanism can be 
simulated by a Coulomb term 

Vo(r) = (32) 

with as the strong coupling constant. This formula must be modified if we consider the 
retardation effects. Indeed, the quark and the antiquark are able to exchange one gluon if 
their separation r is light-like, that is to say if a 2 — r 2 = 0. The probability for a to be 
negative or positive is 

r° l 

p (a < 0) = / daA(a) 2 = - = p(a>0). (33) 

n n 

Consequently, as in the COQM |3j,|4|, we make the following substitution 

V c - V c ^ [8(a + r) + 8{a - r)] , (34) 

where A is an energy scale which is introduced so as to give the correct dimension. This 
parameter is purely phenomenological and could depend on the quark masses i3j]. In this 
paper, we assume that A is a constant. The effective Coulomb potential Vc, treated as a 
perturbation, is then computed with the relation 

/+oo 1 
daA{a) 2 V c — [S(a + r) + S(a - r)] , (35) 

and we obtain a damped effective Coulomb potential, 

^ = -S(^ 1/2 -P(-/"- 2 )^ ^ 



IV. APPROXIMATIONS FOR AM 

The key ingredient to compute the retardation term (j31aj) is the knowledge of p 2 and a 3 , 
respectively given by formulas (j30p and ()A.3|) . As these expressions are complicated, we will 
try to get various simpler ones following the value of the quark mass m. In the following 
we will use two limiting cases for and p 2 quantities: y = corresponds to a vanishing 
angular momentum or to a very high mass, and y = 1 correspond to a very high angular 
momentum or to a very small mass. Both situations are summarized in Table U] 

The solutions of the RSM equations are, in good approximation, very close to the solutions 
of a two-body spinless Salpeter equation with a linear potential 

H ss = 2^p 2 + m 2 + ar (37) 



TABLE I: Values of and p 2 for y = and y = 1. 



y = o 



l 



«3 



p ar 

2 

ap a ap* 
2r 12 6pr 



p irar 

^Xl^ 2 

nap 7r a nap* 



8r 



64 



8/xr 



and with the pure string correction treated as a perturbation |lj|. For massless quarks, this 
correction is only about 6% of the meson mass 16|. So, we will work with the Hamiltonian 
H ss without the string correction, except in Sec. I VI II Actually, this approximation is 
equivalent to consider the RSM at the order y 2 16j|. Within this framework, the extremal 



value of the auxiliary field /i is 



m . 



(38) 



We have then 



and < (p 2 ) < (ji)' 



(Pr) ~ (V) 2 ~ m ' 



(39) 



A. High quark mass 

If we assume that (//) ^> a(r), we can set y ~ 0. In this case, the relations ()31aj) and 
cj) with (fj) ^> a(r) reduces to 



A-^W 1 -®- <40) 



^"-VSv 1 -®- (41) 

It is clear that (p 2 ) is positive. If m is very large, the dynamical effects can be neglected 
and, using (/i) « m, we have 



AM - ^ -»/a^y < 43 > 
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The eigenvalues M are then given with a good accuracy by [17 1 

fa 2 \ 1/?J 

M « 2m + — e n£ , (44) 
\m / 

where e n £ is an eigenvalue of the dimensionless Hamiltonian (<f 2 + |af |). The nonrelativistic 
virial theorem implies that 

2 2 /a 2 \ 1/3 

a(r> = -(M -2m) = - - e n£ . (45) 

3 6 \m J 



So we obtain 



(am) 2 / 3 nr 

/?«» « - — ^ — y — , (46) 



1 /a 2 \ 1/3 /3 




2 Vm / V e«£ 
Approximate values for the quantities e n £ can be found in Ref. 



(47) 



B. Vanishing quark mass 

In this section, we will work with the Hamiltonian H ss for m = 0. With these conditions, 
the relativistic virial theorem gives the following results j^] 

M « (# 5S > = 4(//>, (48) 
o(r> = 2</z). (49) 

Let us first focus on the case £ = 0, for which y = and (p 2 ) ~ (/i) 2 . Thanks to these 
relations, we have for light quarks 

AU = f, (50) 

AM,U = (51) 

Secondly, let us consider the case I ^> 1, for which we can assume that y — 1. We then find 



a / / 7r 16 (p 2 ) 



*|'»' = 4\M 1+ 4-^r)- < 52 > 



ii 



Combining Eqs. and p ty . we find 



(Pi) ~(V) 2 ^ (54) 



Consequently, for massless quarks with high angular momentum, we have 



a fir 64a 2 £(£+l)\ 



a ( ± \ j_(ir 64a 2 £(£+r 



AM '|- ^ -Mo 1 J v * U + ~^r^ » • (56) 



Again, (p 2 ) is positive for £ = and £ ^> 1. 
In good approximation, it appears that 

AM; oc -V. (57) 
M v ; 

The square meson mass composed of light quarks is then given by 

Mf « M 2 + 2M AM,, (58) 

where the term AM 2 is neglected. This shows that the retardation term only causes a global 
shift of the square meson masses and consequently preserves the Regge trajectories, since 
M 2 oc i for large values of I Q, Q- 

V. NUMERICAL RESULTS 

In order to obtain better values for the contribution of the retardation, we will compute 
it at the second order in y. In this case, the quantities a 3 and p 2 become 

(as) - (^ + ^)+i^><2/ 2 >, (59a) 



2 



a{p 2 r ) / a 2 a{pf) a (/ i ) 



(p2) ~ V12 Q(p)(r) + 2(7)) + [do 20(p)(r) I2(r) ) ^ (59b) 
These expressions can be calculated using the relation Q| 

(r) 2 (a(r)/6+(/i)) 2 ' 

In the following, the retardation contribution obtained thanks to Eqs. (|59|) will be called 
"exact" by opposition to the more approximate formulas obtained in Sec. IIVI and it will 
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be simply denoted by AM. The exact contribution is compared to the approximate ones 
in Fig. [T] ((p 2 ) is always positive). In this graph, £ = n = 0, but the qualitative features 
of the curves are the same for other quantum numbers. We can thus determine a validity 
domain for each approximation. In fact, AM; is the best for m < 0.175 GeV (u, d quarks, 
which are commonly denoted n). For 0.175 GeV < m < 4.0 GeV (s, c quarks), AMf, is 
rather good, and for masses larger than 4.0 GeV (b quark), AM/^ works well. As expected, 
the retardation contribution is less important when the quark mass is larger, for which a 
nonrelativistic treatment is more justified. One could ask how the systematic substitution 
(fi 2 ) — > (p) 2 does affect the results. Actually, the values obtained by the two methods differ 
at most by 5%. So, we will maintain our choice, which is to use systematically powers of 

(A*>- 

-0.10 
-0.12 
-0.14 
-0.16 

> 

a -o.i8 
< 

-0.20 
-0.22 
-0.24 
-0.26 

1 2 3 4 5 

m (GeV) 

FIG. 1: Comparison between the exact retardation term and the different approximations of Sec. II VI 
for various quark masses, with a = 0.2 GeV 2 and £ = n = 0. 

We focus now on the light quarks, especially the massless case: The retardation effects are 
indeed expected to be the largest when m = 0, for which the motion is the most relativistic. 
Typical behaviors of AM with I and n are showed in Fig. [21 If we take a = 0.2 GeV 2 , the 
ground state mass M is 1.413 GeV. The contribution of the retardation is —0.205 GeV. 
So, in the worst case, the contribution is about 15% of the non perturbed mass. This result 
justifies a posteriori the perturbative theory we built in Sees. [0] and IIHI Moreover, we 
see that, for a fixed quark mass, the retardation contribution decreases when the different 
quantum numbers, i or n, increase. This means that the key element is not the quark mass 
m, but the its constituent mass (/x) j3] which also increases with these quantum numbers. 
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The more the constituent mass is large, the more the retardation effect is small. 

-0.06 
-0.08 
-0.10 
-0.12 

> 

s -o.i4 

< 

-0.16 
-0.18 
-0.20 
-0.22 

1 2 3 4 5 6 

I or n 

FIG. 2: Exact retardation term versus: £ for n = (filled circles) and n for i = (empty circles), 
with a = 0.2 GeV 2 and m = 0. 

In Sec. IIV[ we showed that the retardation only causes a global shift of the Regge tra- 
jectories for light quarks. Even if this result is only approximate, we can check in Fig. ^ 
that we can have confidence in our formula AM; when m = 0. As a supplementary test in 
this case, we have plotted in Fig. El the square meson masses versus the angular momentum, 
with and without the exact retardation term. We clearly see that the linearity of the Regge 
trajectories is preserved as well as the slope. 

18 
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FIG. 3: Regge trajectories with exact retardation term (circles) and without (triangles), for n = 
and 1, with m = and a = 0.2 GeV 2 . 

Another interesting quantity to study is /3, given by Eq. (|31cjh the size of the relative 
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time part of the wave function (|31b|) and the range of the effective Coulomb potential (|36j). 
Let us see what happens for two extreme cases: the n quark and the b quark, for which we 
take m n = and mb = 4.660 GeV. For n = i = 0, fy, used for the n quark, is given by 
formula (jHOj) . For the b quark, we can use formula (|46jl . For a standard value a = 0.2 GeV 2 , 
we find 

(3hh\ e=n=0 = 0.542 GeV 2 > A| fcn=0 = 0.1 GeV 2 . (60) 

So, the wave function is much more peaked around a = when the mass increases. The 
numerical evaluation with Eqs. (J5T?j) gives 

P\i=n=0 = °- 533 GeV2 > P\t=n=0 = °- 096 GeV2 ' ( 61 ) 

result very close to the approximate one (|5Uj) . 

The parameter (3 also considerably affects the Coulomb potential through Eq. (136)1 . as it 
is shown in Fig.HJ We will see in Sec. IVIII that the effective potential can become very small 
with respect to the retardation term. 

0.0 
-0.2 
? -0.4 



-0.8 
-1.0 

12 3 4 

r(GeV~') 

FIG. 4: Pure Coulomb potential with as = 0.4 (straight line), and effective Coulomb potentials 
for m = and m = 4.660 GeV with A = 1 GeV (dotted lines). 

VI. COMPARISON WITH OTHER APPROACHES 

Since the relative time is introduced in the same way for our model and the COQM, 
both approaches share some common properties. For instance, the relative time part of our 
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wave function Eq. (j31b|) and the counterpart in the COQM have the same form. In this last 
model, for two quarks with the same mass m, we have P, J| 



= ( —J exp ^— — cr 2 ) with /. 



where is a constant related to the interquark potential 

U 



-Kr 2 . 
2 



(62) 



(63) 



Moreover this model predicts also linear Regge trajectories, with a constant square mass 
shift AM^ q M due to the retardation which is given by 



A^coqm = -2V2mK. 



(64) 



Besides this formal similarities, the physical content of these two models are nevertheless 
very different: A comparison is made in Table ITT1 

TABLE II: Comparison between the COQM and our model. 



COQM 



Our model 



Mass operator 
Confinement 
Allowed masses 
Treatment of retardation 



square mass operator 

quadratic 

non zero 

exact 



Time and space decoupling complete and exact 
Shift in energy negative and constant 



ordinary Hamiltonian 

linear 

all 

in perturbation 

partial and approximate 

negative and state dependent 



Relative time wave function gaussian with constant size gaussian with state dependent size 



A more usual method to get relativistic corrections (including retardation) is to consider 
the v 2 lc 2 terms of the Bethe-Salpeter equation or of the Wilson loop formulation of QCD. 

fin n 

This was done for example in Refs. p,|6j. In Ref. j6j, the relativistic correction to the linear 
confinement potential due to the retardation is given by 

a f£(£ + l) 2 



* M = -\^\* +rp ')l- (65) 

As this correction term is obtained by a expansion in v 2 /c 2 , one could expect that the 
best agreement with our method will be obtained for large quark masses. However, these 
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corrections are very different from our term (|43|). In particular, the correction (}4*3*)) decreases 
with the quantum numbers n and I because it decreases with the constituent mass, but the 
contribution (f63|) increases with these quantum numbers. So, our approach does not lead to 
a nonrelativistic limit compatible with previous works. 

Actually, in order to obtain linear Regge trajectories, it is necessary to obtain AM oc 1/M 
(see Eq. (Jo'HjO for light quark systems. But, for heavy quark systems, one could expect that 
AM oc 1/m 2 with m « M/2. Our model clearly misses this transition. A proper treatment 
in a covariant formalism with constraints could probably cure this flaw, but it is out of the 
scope of this paper. 



VII. COMPARISON WITH EXPERIMENTAL DATA 

We saw in Sec. IVII that our formalism had unusual features, compared with already 
known results. So it is important to verify if it can correctly reproduce the experimental 
data. We will make here such an attempt with the nn and bb mesons, in order to check our 
model in different mass domains. The ingredients we pu t in our "realistic" model are: The 
Hannltcnian 63 including the 8tr in g section Qfl, the exact Nation t™ AM, 
the effective Coulomb potential (|36J) treated as a perturbation and the quark self-energy 

Q. 

The string correction has the following form 

AM ■ ^+mM ( ) 

The quark self-energy is due to the color magnetic moment of the quark propagating through 
the vacuum background field, and it has been shown that it brings a contribution to the 
meson mass given in the symmetrical case by 

. , , far)(m/5) . „. 

AM QSE = - J —J±-L±, 67 

with / e [3,4] and 5 e [1.0, 1.3] GeV. The function rj is given, for instance, in Ref. Q], in 
which a more detailed discussion about the quark self-energy and its consequences on the 
meson spectrum can be found. Let us note that 7/(0) = 1. 

The physical parameters we use are given in Table ITTT1 We have tried as much as possible 
to choose standard values: a = 0.192 GeV 2 and as = 0.4 are widely used, and m;, = 
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TABLE III: Our set of physical parameters. 



a = 0.192 GeV 2 a s = 0.4 

m n = 5 = 1.0 GeV 

m b = 4.660 GeV / = 3.0 

A = 1.0 GeV 

4.660 GeV is an acceptable value for the b quark. The parameter / is fixed at 3, which 
is the value obtained by simulations in unquenched lattice QCD calculations [2l| . The 
value 6 = 1.0 GeV is used, but this choice has a very little influence [lfj]. Finally, we fix 
A = 1.0 GeV in order to find, with the above parameters, the nn ground state near the 
center of gravity of the 7r and p mesons, at 612.5 MeV (see below). 

Since our model includes neither the spin (S) nor the isospin (/) of the mesons, the 
experimental data chosen are the spin and isospin averaged masses for the light mesons, 
denoted M av (see Table ITV]) . These are given by [2^ 

E J , J (2/ + 1)(2J+1)M J|J 
" " E^(2/ + 1)(2J+1) ' (68) 

with J = L + S. Mi t j are the different masses of the states with the same orbital angular 
momentum £. For the bb mesons, we present results for the radial excitation of the T (see 

n 

Table EJ. These dat a are taken from Ref. 1231 1 . 

TABLE IV: Comparison between the spin averaged masses M av of some nn family states (see 
Ref. [la] for more details) and the numerically computed masses M (|24|) of our model. The first 
three columns present the different states used to compute the spin averaged masses. The last 
column gives the contribution of the effective Coulomb term. 



Family 


N 2S+1 Lj 


M av (GeV) 


M (GeV) 


(Vc) 


) (MeV) 


P 


l 2S+1 Sj 


0.613 ±0.011 


0.631 




-19 


a 2 (1320) 


l 2S+l Pj 


1.265 ±0.011 


1.235 




-6 


/)(1700) 


l 2S+1 Dj 


1.676 ±0.012 


1.669 




-3 


a 4 (2040) 


l 2S+l Fj 


2.015 ±0.012 


2.022 




-1 


We see 


in Table ITVl that our 


results are in good aj 


^reement with the 


spin averaged 


masses. 
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For all states, the relative error is below 3%. In each case, the influence of the effective 
Coulomb term is very small and could be neglected. Its role in lowering the mass is played 
by the contribution of the retardation. Despite this unusual effect, our approach allows to 
correctly reproduce the spin averaged masses. With a smaller value for the parameter A, the 
contribution of the Coulomb term could be enhanced, but probably to values below those 
obtained in other potential models. 

As it can be seen in Table El the relative error between the data and our result is below 
1% for the mesons T. As expected in this case since j3 is larger, the contribution of the 
Coulomb potential is larger. Heavy quark systems being more sensitive to the very short 
range part of the interaction, better results could probably be obtained by using a running 
coupling constant as{r). But this is out of the scope of this paper. 



TABLE V: Same as in Table EH but for the masses M exp of some bb mesons. The experimental 
error bars are smaller than 10 MeV and are not indicated. 



State 


N 2S + l Lj 


M exp (GeV) 


M (GeV) 


(V C ) (MeV) 


T(15) 


l 3 5i 


9.460 


9.582 


-87 


T(25) 


2 3 Si 


10.023 


9.990 


-52 


T(35) 


3 3 S 1 


10.355 


10.294 


-40 


T(45) 


4 3 5i 


10.580 


10.555 


-33 


T(10865) 


5 3 Si 


10.865 


10.788 


-29 


T(11020) 


6 3 Si 


11.019 


11.002 


-26 



VIII. CONCLUSION 

In this paper, the retardation effects in mesons are taken into account by the introduction 
of a non zero relative time in the rotating string Hamiltonian j?[ Q , following a procedure 
inspired by the covariant oscillator quark model P, 0]. Treated as a perturbation, the 
part of the total Hamiltonian containing the retardation terms is a harmonic oscillator in 
the relative time variable, with an effective reduced mass and an effective restoring force 
both depending on eigenstates of the Hamiltonian independent of the relative time. The 
fundamental state of this oscillator gives the contribution of the retardation to the masses 
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as well as the relative time part of the wave function. The introduction of the retardation 
also affects the Coulomb part of the interaction, which is replaced by an effective damped 
potential. Systems containing two particles with the same mass are only considered, but 
our approach leads to several interesting results. 

In the light quark sector, the contribution of the retardation is not negligible (around 
200 MeV for massless quarks) but it is small enough to justify a perturbative treatment. 
Within this approach, the Coulomb interaction is strongly reduced but the meson masses 
are lowered by the contribution of the retardation. The linearity of the Regge trajectories 
is preserved, which is the most important feature of our model. At last, the relative time 
wave function is a gaussian function centered around zero, which confirms the validity of 
the equal time ansatz in first approximation. 

When the quark mass increases, the contribution of the retardation to the meson masses 
decreases slowly. The relative time wave function becomes more and more peaked around 
zero, as expected. Unfortunately, our model does not lead to a nonrelativistic limit in 
agreement with previous works. This is probably due to the treatment of the relative time 
which is not compatible with a proper elimination of this unphysical degree of freedom O, Q| • 



Nevertheless, when the effective Coulomb potential and the quark self-energy |16( are 
included in our rotating string model with the retardation effects, meson spectra can be 
computed in good agreement with the experimental data. 

This work must be considered as a trial to compute easily the retardation effects in 
mesons. The simplified and approximate approach used here has no firm theoretical basis, 
but it shows that the contribution of these mechanisms to the masses could be non negligible. 
Moreover, the importance of the retardation correction are not controlled by the quark mass 
m but by its constituent-state dependent-mass (y/p 2 + m 2 ) jl^ . 
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APPENDIX: COEFFICIENTS OF THE LAGRANGIAN 



The coefficients used in the Lagrangian (J2J) are defined by 



a 2 = jUi - C(A*i + A*2) + / dp(/3-()u, 

Jo 

a 3 = m (l-C) 2 + // 2 C 2 + / d/3(/3-C) 2 f, 

./o 



a.i 



ci 



c 2 



d/3 t? 2 z/ 



7/ 



d/3 77 ^, 



(A.la) 
(A.lb) 
(A.lc) 
(Aid) 
(A.le) 
(Alf) 



Using the extremal values ()14a|) and ()14b|) of the auxiliary fields rj and u, we can compute 
the following relations 

(A.2a) 
(A.2b) 
(A.2c) 
(A.2d) 





= — arcsins 1 , 






tZ/3 


1 










sVl - 


s 2 + arcsin s 


V 


2ary t 


) 

-2/2 


d/3(J3 


-<> = 


ar 


Vl-s 2 






dp(J3 


-0'* 


ar 


-sVl - 


s 2 + arcsin s 



1 J/1 



-,V2 



In the symmetrical case, we have ii\ — LLi — /•*, mi = m 2 = m, y\ — yi — y, yt = 2y, 
( = 1/2, = 1/2, and fx = fx/2. Equations (|A.1|) are in this case given by 



ar 

ai = 2ix-\ arcsin y, 

y 

a 2 = 0, 
a 3 



arcsin y 



a i 
ci 



a ar f n 

~2~ry i} 1 ^ 1 ~ y2 + arcsin ^) + ft2 ^3 ~y 2 + arcsin y 



■k— vl - y 



2y 2 

c 2 = «^ (-yy/l ~ y 2 4 arc-sin /y | . 



(A.3a) 
(A.3b) 
(A.3c) 

(A.3d) 

(A.3e) 

(A.3f) 
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